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We report on the methods and results of the first dedicated search for gravitational waves emitted during the 
inspiral of compact binaries with spinning component bodies. We analyze 788 hours of data collected during 
the third science run (S3) of the LIGO detectors. We searched for binary systems using a detection template 
family designed specially to capture the effects of the spin-induced precession of the orbital plane. We present 
details of the techniques developed to enable this search for spin-modulated gravitational waves, highlighting 
the differences between this and other recent searches for binaries with non-spinning components. The template 
bank we employed was found to yield high matches with our spin-modulated target waveform for binaries with 
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masses in the asymmetric range 1.0 Mq < mi < 3.0 Mq and 12.0 Mq < m2 < 20.0 AIq which is where 
we would expect the spin of the binary's components to have significant effect. We find that our search of S3 
LIGO data had good sensitivity to binaries in the Milky Way and to a small fraction of binaries in M3 1 and 
M33 with masses in the range 1.0 Mq < mi, m2 < 20.0 Mq. No gravitational wave signals were identified 
during this search. Assuming a binary population with spinning components and Gaussian distribution of masses 
representing a prototypical neutron star - black hole system with mi ~ 1.35Mq and m2 — 5Mq, we calculate 
the 90%— confidence upper limit on the rate of coalescence of these systems to be 15.9 jt^^Ij^q, where Lio is 
10^" times the blue light luminosity of the Sun. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.80.-d 



I. INTRODUCTION 

Currently, there is a world-wide network of kilometer 
scale interferometric gravitational wave detectors that are 
either at or approaching their respective design sensitivi- 
ties. The network includes the US Laser Interferometer 
Gravitational- wave Observatory (LIGO) iU, the British- 
German GEO600 H and the French-Italian Virgo lH. The 
radiation emitted during the inspiral stage of a stellar mass 
compact binary system is thought to be a likely candidate for 
the first direct detection of gravitational waves using these in- 
terferometers iHI^. The initial interferometers will be able 
to search for binary neutron star systems as far as the Virgo 
Cluster, and higher mass binaries which include black holes 
as far as the Coma supercluster The range of merger rates 
consistent with present astrophysical understanding is sum- 
marized in Ref. 101. When binary formation in star clus- 
ters is taken into account with relatively optimistic assump- 
tions, detection rates could be as high as a few events per 
year for initial LIGO id, |^ [l^l- Merger rates derived for 
binary populations in galactic fields consistent with observa- 
tional constraints from the known galactic neutron star - neu- 
tron star systems are highly uncertain but are likely to lie in 
the ranges (at 95% confidence) 0.1 — 15 x lO^^i'^^ ^lo ^"'l 
0.15 - 10 X lO^Vr"^ for black hole - black hole and 
neutron star - black hole binaries respectively iflll IT2I1 . The 
LIGO Scientific Collaboration (LSC) has searched for com- 
pact binaries with non-spinning stellar mass components in 
data collected during the first, second, third and fourth science 
runs (henceforth SI, S2, S3 and S4, respectively )J0, [il, by 
employing optimal matched filtering techniques lfl4ll wherein 
detector data is cross-correlated with a bank of "templates" 
which represent the best current knowledge of the emitted 
waveforms. 

Studies of compact binaries with spinning components ifTsl 
[161 [l7i [isi [19I I20I1 have revealed that general-relativistic dy- 
namical coupling between the spin and orbital angular mo- 
menta (so long as they are not perfectly aligned or anti- 
aligned) will lead to precession of the binary's orbital plane 
which in turn causes a modulation of the observed gravi- 
tational waves' amplitude and phase. The binary's orbital 
angular momentum L, and therefore its orbital plane, and 
the spin angular momenta of the binary's components Si, S2 
will precess about its near constant total angular momentum 
J = L + Si +S2. The gravitational waves observed from 
a binary depend upon the orientation of the binary relative to 
the detector and are strongest along the direction of its orbital 



angular momentum. The amplitude and phase of the gravi- 
tational waves emitted by the binary that will be observed at 
any particular (fixed) location will therefore be modulated by 
the precession of the binary's orbital plane. This precession 
of the orbital plane is nicely illustrated in Ref. fisl l (see Fig. 2 
and the appendix in particular). Figure[T]compares the gravita- 
tional waveforms we would expect to observe from two differ- 
ent binary systems, one consisting of non-spinning bodies and 
the other consisting of spinning bodies. The precession of the 
binary's orbital plane, which is related to the Lense-Thirring 
effect on gyroscopes in curved spacetimes ll2lll . should not be 
confused with the in-plane precession of a binary's periastron 
which occurs in both spinning and non-spinning systems. In 
this search, we consider signals with frequencies of 70 Hz and 
above, corresponding to orbital frequencies of > 35 Hz. We 
are thus sensitive only to the final stages of the binary inspiral. 
By this point, the binary orbit has been circularized due to the 
emission of gravitational waves (see Fig. 5 of Ref. ll22ll ): so 
the precession of periastron degenerates into a secular term in 
the evolution of the phase. 

The statistical distribution of the spins of black holes in 
inspiraling binaries is not well known ll23l l24ll and until re- 
cently the efforts have focused upon developing techniques 
for the detection of binary systems with non-spinning com- 
ponents (for recent reviews see Refs. ll25ll26ll and references 
therein). The presence of amplitude and phase modulations in 
the observed waveforms will reduce our detection efficiency 
when using matched filter templates which do not include spin 
effects |[l3[ll[ll|Ml. These effects are small for low mass 
binaries or binaries with roughly equal component masses, but 
can be significant for high mass or asymmetric systems such 
as neutron star - black hole binaries. 

This paper reports the methods and results of a search for 
gravitational waves emitted during the inspiral of binaries 
consisting of spinning compact objects. This search uses a de- 
tection template family designed to capture the spin-induced 
modulations of the gravitational waveform which could have 
resulted in them being missed by other searches targeted at 
non-spinning systems. This is the first time gravitational wave 
data has been searched for inspiral signals from binary sys- 
tems with spinning component bodies. 

LIGO consists of three detectors located at two sites across 
the US. The LIGO Hanford Observatory (LHO) in Washing- 
ton state consists of two co-located interferometers of arm 
length 4km and 2km and are known as HI and H2 respec- 
tively. The LIGO Livingston Observatory (LLO) in Louisiana 
consists of a single 4km interferometer known as LI . All three 
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FIG. 1 : The gravitational waveforms predicted from the late inspiral 
phase of two different neutron star - black hole systems, one consist- 
ing of non-spinning bodies (upper plot) and the other consisting of 
maximally spinning bodies (lower plot). Both systems are identical 
apart from the spin of their component bodies. Spin-induced preces- 
sion of the binary's orbital plane causes modulation of the gravita- 
tional wave signal and can be clearly seen in the lower plot. 



detectors were operated throughout S3 which spanned 70 days 
(1680 hours) between October 31, 2003 and January 9, 2004. 
The gravitational waves emitted by stellar mass compact bina- 
ries are expected to be at frequencies detectable by LIGO dur- 
ing the final few seconds of the inspiral as well as the merger 
and ringdown stages of their evolution. We analyze S3 LIGO 
data using a detection template family fl9] which efficiently 
captures the amplitude and phase modulations of the signal. 

In Sec. |II]we discuss the evolution of spinning binary sys- 
tems. In Sec. |III]we describe the waveforms that are used to 
model the emission of the target sources we are seeking to 
detect. These target waveforms include modulations to their 
amplitude and phase in order to simulate the effects of spin- 
induced precession of the source. In Sec.|IV]we describe the 
detection template family that we use to search for these target 
waveforms and in Sec. [V] we describe the design and testing 
of the template bank used. In Sec.|VT]we describe the S3 data 
set and summarize the data analysis pipeline. In Sec. lVIII we 
describe various vetoes which were identified as beneficial to 
this search. In Sec. lVIIll we detail results from this search. In 
the absence of a detection we will calculate an upper limit on 
the rate of coalescences using the measured efficiency of our 
search and an estimated population model of the distribution 
of binary systems in the Universe. In Sec.|lX]we perform an 
upper limit calculation based upon the loudest event candidate 
found in our search. Finally, in Sec.|X]we draw conclusions. 
Throughout we shall assume G = c— 1. 



II. EVOLUTION OF SPINNING BINARY SYSTEMS 

We briefly review the current literature regarding the forma- 
tion and evolution of spinning binary systems. The literature 
available focuses mainly on neutron star - black hole (NS-BH) 
binaries (rather than BH-BH binaries). Later we shall show 
that the template bank used in this search is most sensitive to 
binaries with unequal masses such as NS-BH binaries. It is 
likely that the formation of BH-BH and NS-BH (and indeed 
NS-NS) systems are qualitatively similar and that the discus- 
sion here will be relevant to all cases. 

A typical NS-BH evolution would involve two main se- 
quence stars in binary orbit. As it evolves away from the 
main sequence, the more massive star would expand until it 
fills its Roche lobe before transferring mass to its compan- 
ion. The more massive body would eventually undergo core- 
collapse to form a BH, and the system as a whole would be- 
come a high-mass X-ray binary. As the second body expands 
and evolves it would eventually fill its own Roche lobe and 
the binary would then go through a common-envelope phase. 
This common-envelope phase, characterized by unstable mass 
transfer, would be highly dissipative and would probably lead 
to both contraction and circularization of the binary's orbit. 
Accretion of mass can allow the BH to spin-up. It has been ar- 
gued that the common-envelope phase, and associated orbital 
contraction, is essential in the formation of a binary which 
will coalesce within the Hubble time ll23ll . Finally the sec- 
ondary body would undergo core-collapse to form a NS (or if 
massive enough, a BH). Prior to the supernova associated with 
the core-collapse of the secondary body we would expect the 
spin of the BH to be aligned with the binary's orbital angular 
momentum ||23ll. However, the "kick" associated with the su- 
pernova of the secondary body could cause the orbital angular 
momentum of the post-supernova binary to become tilted with 
respect to the orbital angular momentum of the pre-supernova 
binary. Since the BH would have a small cross-section with 
respect to the supernova kick we expect any change to the 
direction of its spin angular momentum to be negligible and 
that the BH spin would be misaligned with respect to the post- 
supernova orbital angular momentum IztIi . The misalignment 
between the spin and orbital angular momentum is expected 
to be preserved until the system becomes detectable to ground 
based interferometers. 

The magnitude of a compact object's spin is dependent 
upon both its spin at formation (i.e., birth-spin) and the spin 
it attains through subsequent accretion episodes. The dimen- 
sionless spin parameter x is given by J/Af ^ where J is the 
total angular momentum of the compact object and M is its 
mass. For a maximally spinning compact object we would 
have X = 1 ^"d for a non-spinning object x = 0. Although 
the estimated birth-spins of NS and BHs are small, simula- 
tions have shown that accretion during a common-envelope 
phase can allow objects to achieve considerable or even near 
maximal spins ||24|] . Due to uncertainties in both the estima- 
tion of birth-spins and modelling of accretion induced spin- 
up predictions of binary's spin population are fairly uncertain. 
The upper bound on a BH's spin is expected to be x ~ 0.998. 
Torque caused by radiation emitted from the accretion disk 
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getting swallowed by the BH counteracts the increase of spin 
caused as the BH accretes mass ll28ll . The upper bound of a 
NS's spin is estimated by calculating the spin which would 
cause it to break up using variety of models for its equation of 
state. The upper limit is estimated to be x ^ 0.7 12^1 . 

Techniques to measure the spin of accreting black holes us- 
ing electromagnetic observations of their accretion disk are 
described in Ref. ifsoll . Using these techniques the spins of a 
black holes in a handful of X-ray binaries have been mea- 
sured, in the case of GRS 1915H-105 the black hole's spin 
was found to be x > 0.98 Isill . Recent observations of spin 
precession through measurement of pulse shapes from binary 
radio pulsars demonstrate misalignment between the orbital 
and spin angular momenta of these systems, see for example 
Ref. IBI . 

For optimal detection of gravitational waves using 
matched-filter techniques we must construct templates that 
represent our best predictions of the signal. These templates 
must model the spin-induced modulations to the waveform's 
amplitude and phase as accurately as possible while still 
resulting in a computationally manageable number of tem- 
plates covering the detectable parameter space. It has been 
shown previously ITtI SEilMl that if spin effects are ne- 
glected when constructing our templates that our detection ef- 
ficiency will decrease and some spinning binary systems will 
be missed. Spin effects are more pronounced when the sys- 
tem's spin angular momentum is larger than its orbital an- 
gular momentum. The Newtonian expression for the mag- 
nitude of the orbital angular momentum of a binary system 
is |Ljv| ~ rjM^^^u}^'^^^ where M = mi + 1112 is the total 
mass of the system, rj = niim2/]\P is its symmetric mass 
ratio and uj is the instantaneous orbital frequency of the sys- 
tem. For given values of M and lu, the orbital angular mo- 
mentum will be largest for binary systems with equal masses, 
TOi = m2- For systems with unequal masses such as NS-BH 
binaries, the orbital angular momentum will be smaller and 
the spin angular momentum will play a more significant role 
in the system's evolution. It will therefore be more suscepti- 
ble to the effects of spin than equal mass systems (see Fig. 3 
of Ref. lfl6ll ). For schemes that fail to take into account spin 
effects, detection efficiency will be worse for binaries with 
i) unequal mass components, ii) components with large spin 
magnitude and iii) when there is significant misalignment be- 
tween the spins and the orbital angular momentum. 



III. TARGET WAVEFORMS 

In this section we describe the fiducial target waveforms 
used to represent the gravitational wave signals expected from 
binary systems of spinning compact objects. We adopt the 
post-Newtonian (PN) equations given out in Ref. 111911 and 
based upon Refs. 1 15, 16, 33, 34, 35, 36,l3l[3l (see Ref. ^ 
for a complete list of references to all original derivations), 
which model the inspiral of the binary in the adiabatic limit. 
In this limit the binary's components follow a sequence of 
shrinking instantaneously-circular orbits in a precessing or- 
bital plane. 



The instantaneous orbital frequency lo evolves according to 
Eq. (1) of Ref. Hit], which has the structure 

-^ = ^;,(a;(i),L^W-Si,2(t),Si(i).S2(t);A/,r7,xi,2), 

(1) 

with the total mass of the system M, the symmetric mass ratio 
?7, the magnitudes of the binary's dimensionless spin parame- 
ters X1.2, the direction of the Newtonian angular momentum 
L7v(0 (oc r X V, perpendicular to the bodies velocity and 
the vector joining them), and the directions of the two spins 
Si, 2(0- Orbital PN effects are included up to 3.5PN order, 
while spin effects are included up to 2PN order. 

The two spins and the orbital angular momentum evolve ac- 
cording to standard general-relativistic precession equations, 
which are truncated consistently at the relevant PN order, and 
which have the structure 

Si = F^^(w,Lw,S2;M,?7,x2) X Si, 

^2 = Fjcj,LAr,Si;M,r;,xi) X S2, (2) 

Ltv = (tJ,LAr-Si,LAr-S2,Si,S2;M,77,Xl,X2) X Lat 

(see Eqs. (2), (3), and (9) of Ref. |l9ll ). 

The gravitational strain perturbation /i*-' is computed from 
the leading-order mass-quadrupole term specialized to circu- 
lar orbits, following Finn and Chernoff ||39ll (see also Sec. II 
C of Ref. ill). Since Finn and Chernoff use a fixed source 
coordinate system, the twice-differentiated mass-quadrupole 
tensor Q*-' is a function of the orbital phase J ujdt and of 
Ltv {t)- The response of a ground-based interferometric detec- 
tor is obtained by projecting Q^P onto a combination of unit 
vectors along the interferometer arms, which introduces a de- 
pendence on five angles that describe the relative direction (0 
and Lp, which subsumes the initial orbital phase of the binary) 
and orientation (6, (f>, and 1/)) between the detector and the 
Finn-Chernoff source frame. 

Equations ^ and (|2]i are integrated numerically in the 
time domain until the minimum of the PN orbital energy 

■E^3PJv(t^jv,Si,S2,A/, ?/, Xi,X2) (see Eqs. (11) and (12) 
of Ref. 11911) is reached or until cu becomes negative. No at- 
tempt is made to describe the waveform beyond this stopping 
point, where it is assumed that the adiabatic approximation 
must break down. Altogether, the waveforms are functions of 
four mass and spin constants (M, 77, xi, and X2), of six angles 
describing the orientations of Lyv, Si, and S2 at a fiducial 
time and frequency, the five direction and orientation angles 
and the distance of the detector from the source. We note that 
the angles 8 and Lp are degenerate with the angles given im- 
plicitly when we define Lat. In this analysis we assume that 
the binary's orbits have become circularized (see brief discus- 
sion in Secini and that the orbital eccentricity is zero. Given 
this assumption we are able to describe the binary using 15 
independent parameters. 
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IV. DETECTION TEMPLATE FAMILY 

As discussed in Sec. [Ill when the binary components carry 
significant spins which are not aligned with the orbital angular 
momentum, spin-orbit and spin-spin couplings can induce a 
strong precession of the orbital plane, thus causing substantial 
modulation of the gravitational waves' amplitude and phase 
(see Fig. [T]). Detection-efficient search templates must ac- 
count for these effects of spin. A straightforward parametriza- 
tion of search templates by the physical parameters that af- 
fect precession results in very large template banks, which is 
computationally prohibitive. It is then necessary to reduce the 
number of waveform parameters while still efficiently cover- 
ing the parameter space of target waveforms. 

We shall denote by "detection template family" (DTF) a 
family of signals that captures the essential features of the 
true waveforms, but depend on a smaller number of parame- 
ters, either physical or phenomenological. At their best, DTFs 
can reduce computational requirements while achieving es- 
sentially the same detection performance as true templates. 
However, DTFs can include non-physical signal shapes that 
may increase the number of noise-induced triggers, affecting 
the upper-Umit studies. Moreover, DTFs are also less ade- 
quate for parameter estimation, because the mapping between 
template and binary parameters is not one-to-one. 

In recent years several DTFs for precessin g co mpact bina- 
ries have been proposed H [H [11 [H SMSa A 
DTF based on the so-called Apostolatos ansatz OlSl llTIl for 
the evolution of precession frequency was thoroughly investi- 
gated in Refs. lEol l40ll . It was found that the computational 
requirements of the Apostolatos-type families are very high, 
and its signal-matching performances are not very satisfac- 
tory. An improved version using spiky templates was then 
proposed in Ref. iIztIi . 

After analyzing the physics of spinning-binary precession 
and waveform generation, the authors of Ref. [12J showed 
that the modulational effects can be isolated in the evolution 
of the two gravitational wave polarizations (i.e., h+ and /ix), 
which combined with the detector's antenna patterns yield its 
response. As a result, the detector's response can be written 
as the product of a carrier signal and a complex modulation 
factor, which can be handled using an extension of the Apos- 
tolatos ansatz. More explicitly, the modulated DTF in the fre- 
quency domain proposed in Ref. lfl9ll reads: 



e2-/*<'0(/,„t - /) (for / > 0) 



(3) 



with h{f) = h*{—f) for / < 0. The coefficients aj in Eq. (O 
are six real coefficients encoding the global phase, the strength 
of the amplitude modulation, its relative phase with respect to 
the leading order amplitude, and the internal (complex) phase 
of the modulations. The coefficient is the time of arrival 
and 9{...) is the Heaviside step function which is zero for all 
frequencies / > /cut- We use the parameter /cut to termi- 
nate the template waveform once we believe it is no longer 



an accurate representation of the true gravitational waveform 
(generally due to deviation away from the adiabatic approxi- 
mation). 

In Eq. (O the functions hj{f) = Aj{f)e''f'^"'^f^ are the 
basis-templates where Aj (/) are the real amplitude functions; 



(4) 
(5) 
(6) 



where B — i3f~^^^ and (3 is related to the frequency of pre- 
cession ll4lll and is used to capture the spin-induced mod- 
ulation of the waveform. The function ?/'nm(/) represents 
the phase of the non-modulated carrier signal; it depends on 
the masses and spins of the binary's components and it can 
be computed in post-Newtonian (PN) theory. Here, as in 
Ref. 111911 . we express V'nm in terms of only two phenomeno- 



Mf) 


= r 


-7/6 




Mf) 


= r 


-7/6 


cos(S), 




= r 


-7/6 


sin(S), 



logical parameters i/jq and ^3 ||6C 

^nm(/) = /-'/'(V^o 



^3/) . 



(7) 



In the case of single-spin binaries (i.e., only one of the bodies 
has spin), it is possible to (analytically) relate the three phe- 
nomenological parameters V'Oj^ and /3 with the four physical 
parameters M, 77, ki and xi 114111 . The physical parameter k 
is the cosine of the angle between the direction of the (total) 
spin and the orbital angular momentum and in this case would 
be Ki = Ljv- Si. However, for double-spin binaries — which 
is the case investigated in this paper — the mapping is not an- 
alytical and the number of physical parameters is greater than 
four, resulting in an intractably large template bank. Within 
the spirit of DTF and, as a first step in implementing search 
templates for spinning, precessing binaries, we proceed here 
with the three phenomenological parameters ^q, V'3 and /3. 

The DTF described by Eq. (O generalizes the Apostolatos 
ansatz in two ways: it allows a complex phase offset between 
i) the leading order /^^/® amplitude term (Eq.|4|i and the sinu- 
soidal amplitude terms (Eqs.|5]and|6]l and ii) between the co- 
sine and sine modulation terms. Quite interestingly, as shown 
in Ref. iiUl . by an appropriate choice of the phenomenologi- 
cal coefficients ai - e, the DTF also has the ability to generate 
higher harmonics which arise in the target signal discussed in 
Sec.|lII] Those higher harmonics are caused by oscillations in 
the components of the gravitational wave polarization tensor 
and not directly by the precession of the orbital angular mo- 
mentum and spins, and should be reproduced by the search 
templates in order not to lose efficiency. 

Henceforth, we will treat i/jq, V'3 ™d [3 as intrinsic param- 
eters and the ai...Q and to as extrinsic parameters. Intrinsic 
parameters describe the source itself (e.g., masses, spins). To 
maximize the SNR with respect to the intrinsic parameters 
we must construct templates corresponding to different val- 
ues of the intrinsic parameters and measure the SNR obtained 
by each of these templates with our detector data. On the 
other hand, extrinsic parameters describe the observer's rela- 
tion to the source (e.g., distance of the source from the ob- 
server, the amplitude and time of arrival of the gravitational 
wave at the observer). Maximization of the SNR with respect 
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to extrinsic parameters can be performed automatically (e.g., 
measurement of a signal's time of arrival using an FFT) and is 
computationally cheaper than maximization of the SNR with 
respect to the intrinsic parameters. 

In practice we set /cut to the frequency of the gravitational 
wave emission at the last stable orbit (LSO) which we estimate 
using 



/c 



/lso 



3/2 



(8) 



where tlso = 671/ is the separation of the binary's com- 
ponents and the total mass M is estimated from tpo and ?/'3 
using approximate relationships between phenomenological 
and physical parameters we introduce in the next section, see 

Eqs. (fnira . 

To assess whether a stretch of detector data contains a grav- 
itational wave signal we calculate the signal-to-noise ratio 
(SNR) which is the cross-correlation of our templates with 
the data. The full process of deciding whether a detection has 
been made is described in S ec. l VII of this paper and more fully 
in the companion papers lf7l l26ll . We can simplify the calcula- 
tion of SNR by orthonormalization of the amplitude functions 
Ak- We obtain the orthonormalized amplitude functions, de- 
noted Ak, using the Gram-Schmidt procedure which leads to 
the transformations: 



Ai Ai 



Ai 



Ao 



A. 



A2 = 



A2 - (A2,Ai)Ai 



\\A2-(A2,AAAi\\^/^ 



(9) 



^3 = 



^3 - (A3,Ai)Ai - (A3,A2)A: 



11^3 



A3,Ai)Ai - (^3,^2)^2 



where we use | |a| | to represent the inner product of a function 
with itself: ||a|| = {a, a). Throughout we will use the real- 
valued inner product: 



(a, 6) = 45ft / df 



ra *{mf) 

Shif) 



(10) 



where Sh (./) is an estimate of the noise power spectral density 
of the data. The final form of the orthonormalized amplitude 
functions are very long and for that reason not reproduced 
here. The DTF in terms of the orthonormalized amplitude 
functions has the exact same form as that shown in Eq. (O 
with h, hj and aj replaced by h, hj and aj respectively. De- 
manding templates normalized so that {h, h) = (h, = 1 
leads to the constraint = 1- Having defined the or- 

thonormalized amplitude functions Ak we can calculate the 
SNR, p: 



p = max (x, h(to, aj)) 



to 



\ 



^(x,/lj(to)) , (11) 



i=i 



where x is the detector data and the orthonormalized basis- 
templates are given by 

hj = Ajif)e''f'™'^f^ for j = 1, 2, 3 and 

/ij = iA,-3(/)e"^™^-''^ for j = 4, 5, 6. (12) 

Note that we do not explicitly need to calculate ai...Q in or- 
der to calculate the SNR but that they can be found simply if 

required: aj = (^x, hj{to)'j /p. 

For Gaussian white noise, p^ will, in general, have a 
distribution with 6 degrees of freedom. In the case where the 
spin parameter = we find that A2 and ^3 both vanish 
and that p^ is described by a distribution with 2 degrees of 
freedom. To reflect the increased freedom we choose a higher 
SNR threshold, = 12 when /3 ^ and a lower value of 
« 11.2 when (3 — Q. These values were chosen to give 
approximately the same number of triggers when analyzing 
Gaussian white noise and to ensure that the number of triggers 
produced during the real search was manageable. 

V. TEMPLATE BANK 

Since we will not know the parameters describing an inci- 
dent gravitational waveform a priori, we must filter our de- 
tector data with a set of templates known as a template bank. 
Neglecting the effects of noise, we would expect that the tem- 
plate yielding the largest SNR to be the best representation 
of an incoming signal. Due to the discrete nature of the tem- 
plate bank (it must be discrete since it can only contain a finite 
number of templates) we will lose SNR due to mismatches be- 
tween the intrinsic parameters of any gravitational wave signal 
and the best template. By placing templates with an appro- 
priate density, we can limit the maximum mismatch between 
signal and template intrinsic parameters and hence limit the 
loss of SNR caused by the discreteness of the bank. The spac- 
ing of templates in the intrinsic parameter space required to 
limit this mismatch can be found using the metric on the sig- 
nal manifold ll42[ I43I1 . In this section we describe the cal- 
culation of the metric, the template placement algorithm and 
comparisons with other banks before discussing the testing of 
the bank using software-injected simulated signals. 

A. Metric calculation 

In this search we use a simple metric based on the strong 
modulation approximation described below. The rationale is 
that systems with waveforms only weakly modulated by spin- 
induced precession should be detectable with high efficiency 
by a non-spinning binary search, e.g., |01. Thus we concen- 
trate on designing a bank that will capture systems whose 
waveforms will be strongly modulated. The metric calculation 
and template placement (or tiling) algorithms become much 
simpler in the strong modulation limit. More recently, more 
precise treatments of the full metric on the DTF parameter 
space have become available and work is in progress 

to incorporate them into future searches. 
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In the strong modulation approximation, the orbital plane 
is assumed to precess many times as the gravitational wave 
sweeps through the LIGO band of good sensitivity. Also the 
opening angle between the orbital and spin angular momen- 
tum is assumed to be large, corresponding to large ampli- 
tude modulations of the signal. Mathematically this corre- 
sponds to the statement that the precession phase B sweeps 
through many times 27r and thus that the basis-templates hj 
are nearly orthonormal (without need for the Gram-Schmidt 
procedure). Below we shall see that this assumption places a 
condition on the precession parameter /?, which for the initial 
LIGO design noise power spectral density ll45ll corresponds to 
P > 200 Hz2/3. 

We can relate this condition for validity of the strong mod- 
ulation approximation to the astrophysical parameters of sys- 
tem. Naively we can put the phenomenological parameters in 
terms of astrophysical parameters using: 



-00 



3 r / , M-5/3 {mi+m2f 
— [7r(mi +777,2)] ^'^ — zr-Z 



777 1 7772 



, , M 9 (7771 + 7779)^ 



777 1 7772 



(13) 



(14) 



Mr 



Ami J 7772 V"^l + "^2 



2/3 



(15) 



acknowledging that, in reality, the true signal manifold and 
phenomenological template manifold do not map this simply. 
The equations for t/jq ™d 7/^3 can be found by considering 
the expansion for the gravitational wave phase_!/'(/) given in 
terms of masses (e.g., Eqs. (3.3) and (3.4) of i4g]) and equat- 
ing the dominant terms of this expansion to those with the 
same frequency exponent in the expansion for gravitational 
wave phase given in terms of t/'o.i,... in lfl9ll47ll . The effects 
of spin are neglected in these approximations of V'o and 7/13. 
The equation for j3 arises by recognizing that (3 is related to 
the evolution of the rate of precession, see Eq. (45) of ifTsIl and 
ll4lll for further discussion ll6lll . 

The constraint for validity of the strong modulation approx- 
imation is that the mass ratio must satisfy 1112/ rrii > 2. Also, 
we specify that the total mass be less than some value (here 
^15 Mq) so that the waveforms do not begin far enough into 
the non-linear region to require extra phenomenological pa- 
rameters. Thus the parameter space region of such a search 
may be expressed solely in terms of the range of masses for 
the lower-mass body. In this search the range used for 7T7i 
was 1.0 Mq < nil < 3.0 Mq, a likely range of masses for 
neutron stars, and 6.0 Mq < m2 < 12.0 Mq was used as 
the range for the more massive body. Thus, astrophysically 
this search is directed at NS-BH systems or BH-BH systems 
with unequal masses. These mass ranges are converted into 
ranges of ipo and 7/^3 using Eqs. (T3[ and ( fT4l i which define the 
region of parameter space we populate with templates. Due to 
the inexact nature of these equations we know that the range 
of masses for which the template bank obtains its highest 
matches may differ from the range of masses we use to specify 
the region of {tpo, tp^). We will show in Sec. lVCI that a tem- 
plate bank generated using the range of masses just specified 
yields high matches (greater than 0.9) for binaries with physi- 



cal masses in the asymmetric range 1.0 Mq < 7711 < 3.0 Mq 
and 12.0 Mq < 7772 < 20.0 Mq. We shall also show that 
this search is efficient for non-spinning systems as well as for 
spinning ones. 

We derive the metric components in the manner of 
Ref. ll42ll . Starting from the detection statistic (the square 
of Eq. ([TT])). let us take our data x to have the form of a tem- 
plate with parameters slightly perturbed from those of the tem- 
plate h we filter it with: 

= (ai + M2)e*('^™+'*'^™) 

+ (a3 + iai) cos{B + dS)e'W'N"+'''^™) 

+ (a5 + iae) sm{B + rfe)e*('^NM+rf^NM) (jg) 

Note that only the intrinsic parameters are perturbed, as the 
maximization takes care of the extrinsic parameters. Expand- 
ing to second order in the perturbation, we have 



+ (Q!3 + 7^4 ) 

+ (q;5 + iaa) 



dB' 
dB^ 



(ai 
h2- 

/73H 



h ia2)hi 
dBhs 
dBh2 



(17) 



Under the approximation that the hj are orthonormal, we get 



(x, hi) 
{x, hi) 
(x,h2) 

{x,h^) 

(.T, h^) 

{x, he) 



ai 



a2 



1 - (d^lu) 



-F {di, 



2 1 



- a2-F(d7/;NM), 
+ aiF{d%j}^M), 



1 - \f {d^U) - \f {dB' 



0:4 



-a4i^(d7/'NM) + a5F{dB) - ae^C^V'NM dB) 

22 



1 - If {d^b^) - If {dB' 



rj V " T- IN m / rj 

aeF{dB) + a^Fidip^M dB), 



"5 



1 



1 



-asFidB) + a4F(d7ANM dB), 



1 



^F{di:lu)-iF{dB') 



2 ^-'^^^"^^ 2 
aiF(dB) — a^Fidip^M dB), 



+ a3F((i-0NM) 

- a6F((iV'NM) 

+ a^F{dipi^M) 
(18) 



where F is a functional (originally defined in Ref. i42|| as J) 
given by 



/„ax//0 



-7/3 



F{a) = — / dx 
and the noise moment / is itself defined as 



a{x) 



/max//0 ^-g/3 



(19) 



(20) 
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where /min and /max define the range of frequencies we inte- 
grate over. In S3 we used a lower cutoff frequency of 70 Hz, 
chosen to exclude lower frequencies for which the detector's 
power spectral density was significantly non-stationary, and 
an upper frequency corresponding to the Nyquist frequency, 
in this case 1024 Hz. Inserting the relations from Eq. (fTsT l 
into Eq. (fTTT i and keeping up to second order perturbations, 
we obtain 

6 6 

- Y.a][F{dB^)~F{dBf] 

j=3 



hyM 



2 (aattg - a4a5) 



[F(rf^'NM dB) - F{d^^M)F{dB)] 



(21) 



To finish computing the perturbed we must maximize 
Eq. ( l2n i over the coalescence time and aj (subject to the con- 
straint X]j=i '^^ ^ ^ since we are dealing with normalized 
waveforms). Maximization over aj is performed straightfor- 
wardly using Lagrange multipliers. We find ai ~ a2 ^ 0, 
Q!3 = — ag, and a4 = as, which leads to 

max {x, hjf = 1 - Fiddly,) + F(dV'NM)' - ^(dS^) 

+F{dBf + F(dV'NM dB) - F{di'j^M)F{dB). (22) 

We incorporate the time-dependence of into the template's 
phasing and expand the phase functions in terms of the phe- 
nomenological parameters and coalescence time tc 





-5/3 + rf^3 /-2/3 ^ 27r/dt„ 


(23) 


dB = dpf- 


2/3 


(24) 


Using the definition of the metric ll42ll to write 






1 - 2,9,fcdA'^dA^ 


(25) 


we obtain the metric components 






4^2 (J^ _ J2) ^ 




25tciAo = 


27r ( Jg - J4 J12) , 




25tciA3 = 


27r ( Jg — J4 Jg) , 




25tc/3 = 


(-7r/2) ( Je - J4 Jg) , 




2(7i/'oi/'o ^ 


Jl7 - Jl2: 




2fl'l/'0'/'3 ^ 


Ji4 — JgJi2, 




25^-0/3 = 


(-l/2)(Ji4- JgJi2), 






Jll - Jg, 




25^3/3 = 


(-1/2) (Jll - J|) , 




2ff/3/3 = 


Jll — Jg 


(26) 



before projecting out the coalescence time tc- Here we have 
used Jg to represent the normalized noise moments given 



J, = I,/l7 



(27) 



where the noise moment / was defined in Eq. ( |20l i. These mo- 
ments give us a way of checking when the strong modulation 
approximation is valid. 




100 200 300 400 

P (Hz^^^) 



500 



FIG. 2: Moment functions CriP) (solid line) and Sr{l3) (dashed 
line) for the initial LIGO design noise power spectral density. For 
values of 13 > 200 Kz^^^ we see that these moments become small 
and can be neglected - this is what we call the strong modulation 
approximation. 



If we had not made the strong modulation approximation, 
we would also need the functions 



/ df 

Jo 




\osB{f)/h, 


(28) 


/ df 
Jo 




'\mB{f)/lT, 


(29) 



which we call the cosine and sine moment functions. The in- 
ner products of the basis templates hj with each other (prior 
to the Gram-Schmidt procedure) are proportional to these mo- 
ment functions, and thus the strong modulation approximation 
corresponds to assuming that Cj and 87 are small compared 
to unity. For the initial LIGO design noise power spectral den- 
sity curve ll49ll the moment functions are plotted in Fig.|2] We 
see that the strong modulation approximation should hold (to 
about the 10% level) for (3 > 200 Hz^/s. See also Fig. 15 of 
Ref. ll4lll . discussed more below, which shows approximately 
the same behavior 



After projecting the coalescence time out of Eq. ( |26l ) and 
dropping i/;/3 cross terms (which simplifies the template place- 
ment and changes the volume per template by less than 3%), 
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we obtain 

'^9i>oi>o = '^IT — Jl2 ^ ^ JiJu)^ I {Jl - Jf) , 

2ff/3/3 = Jii - J| - ( Je - J4 J9)V4 (Jl - J|) • (30) 

B. Template placement algorithm 

We set the density of our template bank in terms of the min- 
imal match [M M), defined to be the lowest match that can be 
obtained between a signal and the nearest template ||42||. A 
template bank designed to have minimal match MM — 0.95 
would therefore suffer no more than a 1 — Af M ~ 5% loss 
in SNR due to mismatch between the parameters of a signal 
and the best possible template in the bank (assuming that the 
signal and templates are from the same family). 

The metric components shown in Eq. (|30] | are constant in 
the strong modulation approximation, which enables us to 
use a simple template placement algorithm. We use a body- 
centred cubic (BCC) lattice which is the most efficient tem- 
plate placement in three dimensions. We first diagonalize the 
metric, which leaves the (3 parameter unchanged but gives us 
new "horizontal" parameters ?/'q and i/jg. Starting on the plane 
/3 = 0, we draw a box in the primed coordinates which en- 
closes the part of that plane to be searched. Beginning at one 
corner of this box, we step in the primed "horizontal" coor- 
dinates by amounts (4/3) a/2 (1 - MM)/E, where E is the 
corresponding eigenvalue of the metric, i.e., g,p'^,p'^ or 
At each point we transform to the mass parameters using 
Eqs. ( fTlT l and ( fT4l ) and check if we are in the targeted region of 
physical mass space. If the point is within that region, we add 
a template to the list. Once a plane of constant {3 is filled, we 
move "up" a distance in /? equal to (2 /3) ^2(1 — MM)/g0(3, 
and lay a "horizontal" grid which is staggered half a cell (in 
both primed directions) from the previous one. Thus a BCC 
lattice is formed. 

Such a simple template placement algorithm is susceptible 
to the "ragged edges" problem. That is, there will be some 
areas near the edge of the targeted region of parameter space 
that will match the nearest template at a level less than AIM. 
The problem appears in other template placement algorithms 
such as those of Refs. Izsi l43ll . and sometimes is addressed 
in a complicated way. Our solution is simple and practical. 
In stepping around the {ipQ, 1/^3) plane, we check to see if we 
have crossed the edge of the targeted region. If we find our- 
selves at a point outside of the targeted region, we check to 
see whether the point halfway between the current position 
and the previously laid template is itself within the targeted 
region. If so, we add a template there. Although the edges of 
the targeted region are curved, the radius of curvature is many 
template spacings meaning that we can treat the edges as fairly 



straight. This simple method solves the ragged edges problem 
while resulting in a small number of additional templates. 

As mentioned earlier, we choose /cut, the frequency at 
which we end our template, to be the frequency of gravita- 
tional wave emission at the last stable orbit. However, we 
compute metric components by effectively taking /cut to in- 
finity, which gains us simplicity at the cost of a small over- 
coverage. 

We can compare the simplified template bank used here to 
those proposed in the literature, particularly in Refs. iiHliill . 
Although neither of those articles actually constructs a tem- 
plate bank or gives explicit metric components, we can find a 
point of comparison. Figure 15 of Ref. ll4lll plots the coordi- 
nate volume per template as a function of /3, assuming a sim- 
ple cubic lattice with MM = 0.97 and an analytical approxi- 
mation to the initial LIGO noise curve. In the high-/3 (strong 
modulation) limit, their volume tends to ~ 5 x 10^ Hz'^. For 
the same M M, lattice, and noise curve, our volume per tem- 
plate is ~ 6.4x10® Hz"^. Thus, our grid is slightly sparser than 
that of Ref. ll4lll . Most of the difference is because they define 
their final metric (on the space of intrinsic parameters only) in 
terms of a "minimax" overlap, which is more restrictive than 
the metric described here. The issue is that the spacing on the 
intrinsic parameter space in general depends on the extrinsic 
parameters, and there are multiple ways to remove this depen- 
dence. The minimax criterion of Ref. iiHl assumes the worst 
case (in terms of extrinsic parameters) or tightest spacing for 
each point in parameter space, and thus is tighter (lower tem- 
plate volume) than it needs to be. Spin-induced precession of 
the orbital plane will cause sidebands either side of the carrier 
frequency. The metric we describe is constructed implicitly 
assuming that there is always non-zero power at the carrier 
frequency and both precession sidebands, which eliminates a 
set of measure zero of worst-case points in the extrinsic pa- 
rameter space. The template bank tests described below ver- 
ify that the loss of efficiency due to neglecting the worst-case 
extrinsic parameters is no more than a few percent. 

For the real S3 noise spectra which were used to construct 
the template banks in this search, template numbers were typ- 
ically 2 — 6 X 10"^ in HI and LI when prescribing a minimal 
match of 0.95. The number of templates was larger in H2 
compared to the other detectors and also increased with time 
to ~ 1.6 X 10* towards the end of S3 due to a flattening of the 
noise power spectrum in H2. Although a minimal match of 
0.95 was prescribed the effective minimal match of the tem- 
plate banks generated was reduced to 0.93 due to a small 
calculation error Figure [3] shows a template bank generated 
using 2048 seconds of HI data and with a prescribed minimal 
match of 0.95. 



C. Testing the template bank 

The template bank was tested using a series of simulated 
signals constructed using the equations of the target wave- 
forms described in Sec. |III] We considered a variety of spin 
configurations including systems where neither, one or both 
bodies were spinning. We also considered masses outside 



11 




BOO 
CO- 400 

200 




600 
CO. 400 

200 



-5000 -4000 -3000 -2000 -1000 



FIG. 3: A template bank generated with minimal match =0.95 us- 
ing 2048 seconds of HI data taken during S3. The crosses show 
the positions of individual templates in the {tpo , tps , f3) parameter 
space. For each template a value for the cutoff frequency /cut is es- 
timated using Eq. [8] This bank requires a 3-dimensional template 
placement scheme in order to place templates in the {ipo, ^ps, 13) pa- 
rameter space. Previous searches for non-spinning systems have used 
2-dimensional placement schemes. 



the range we expected the template bank to have good cov- 
erage in order to fully evaluate the range of masses for which 
it could be used. For each spin configuration we created a 
series of signals corresponding to every mass combination: 
1.0 Mq < mi,m2 < 20.0 Mq. Using the initial LIGO de- 
sign sensitivity we then measured the best match that could 
be obtained for every signal using our template bank. Fig- 
ure |4] shows a sample of the results from the tests of the tem- 
plate bank. As expected we found that our template bank 
achieved the highest matches for non-spinning (and there- 
fore non-precessing) binaries. Performance degrades as spin- 
precessional effects become more pronounced i.e., when both 
bodies are spinning maximally with spins misaligned from 
the orbital angular momenta. The template achieved matches 
> 0.9 for a mass range 1.0 Mq < mi < 3.0 A/0 and 
12.0 Mq < 7712 < 20.0 Mq (and equivalent systems with mi 
and m2 swapped). The detection template family (described 
in Sec. lIVI i is capable of obtaining high matches for compara- 
ble mass systems, the lower matches obtained for comparable 
mass systems are a result of targeting our template bank on 
asymmetric mass ratio systems (which are more susceptible 
to spin effects and conform to the strong modulation approxi- 
mation). 

Matches below the specified minimal match of 0.95 in the 
bank's region of good coverage are a consequence of (small) 
differences between the DTF and the target waveforms mean- 
ing that the DTF cannof perfectly match the target waveforms. 
The fitting factor (FF) measures the reduction of SNR due to 
differences between the DTF and the target waveform ITtIi 
(and should not be confused with the minimal match which 
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FIG. 4: Plots showing the best match achieved by filtering a series 
of simulated signals through the template bank described in this sec- 
tion. The values on the x and y axes correspond to the component 
masses of the binary source to which the simulated signal corre- 
sponds. The shade of grey in the plots shows the best match achieved 
for a given simulated signal. The four subplots correspond to four 
different spin-configurations of the binary source. The top-left sub- 
plot shows results for a non-spinning binary system. The top-right 
subplot shows results for a system consisting of one non-spinning 
object and one maximally spinning object with its spin slightly mis- 
aligned with the orbital angular momentum. We would expect this 
system to precess. The bottow two subplots show results for two 
generic precessing systems consisting of two maximally spinning 
bodies with spins and orbital angular momentum all misaligned from 
each other. We see that the region of the mass plane for which we ob- 
tain matches > 0.9 is largest for the non-spinning system and tends 
to be concentrated in the asymmetric mass region loosely bounded 
by 1.0 Mq < mi < 3.0 Mq and 12.0 Mq < 7712 < 20.0 Mq. 



measures the loss of SNR due to discreteness of the template 
bank ll42ll ). The DTF performance is evaluated and its fitting 
factor is measured in Sec. VI of Ref. 11191], fo r NS-BH systems 
an average FF of « 0.93 was measured 116211 . 



VI. SEARCH PIPELINE 

The pipeline used for this search is the same as used in the 
other S3 searches for binary inspirals |01 and is described fully 
in a set of companion papers ll26ll50ll . This pipeline has been 
significantly updated since the S2 analysis and a brief sum- 
mary is now given. 

In Sec. I VI Al we discuss the S3 data set. In Sec. lVIBl we de- 
scribe how we decide whether triggers measured in different 
detectors could be associated with the same gravitational wave 
event. In Sec. lVICl we introduce the statistic which we use to 
assign SNRs to the events found in coincidence between two 
or more detectors. In Sec. |VIDl we describe how we estimate 
the expected rate of accidental coincidences. 
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A. Data sample 

To begin with we construct a list of times for which two or 
more of the detectors are operating nominally, in what is re- 
ferred to as science mode. By demanding that a gravitational 
wave be detected in coincidence between two or more detec- 
tors we simultaneously decrease the probability of inferring 
a detection when no true signal was present (a false alarm) 
and improve the confidence we have in a detection of a true 
signal. Data collected by the LHO detectors was only ana- 
lyzed when both detectors were in science mode. This was 
due to concerns that since both of these detectors share the 
same vacuum system, the laser beam of a detector in anything 
but science mode might interfere with the other detector. 

We denote periods of time when all three detectors are in 
science mode as H1-H2-L1 times and periods when only the 
Hanford detectors are on as H1-H2 times. A coincident trigger 
consisting of a trigger in the HI detector and the LI detector 
wiU be referred to as an Hl-Ll coincident trigger and simi- 
larly for other combinations of detectors. 

In this search we analyze 184 hours of H1-H2-L1 data and 
604 hours of H1-H2 data (see Table|I]i. During these times we 
construct template banks for each detector and subsequently 
produce a list of triggers whose SNR exceeded our threshold. 

Around 9% of the data is specified as playground data and 
is used to tune the various parameters (e.g., SNR thresholds 
and coincidence windows) used in the full search. Playground 
data is not included in the upper limit calculation but is still 
searched for possible detections. We also construct lists of 
veto times during which the data we analyze had poor data 
quality due to short stretches of instrumental or environmen- 
tal noise ifsol Isill . All coincident data is analyzed but gravi- 
tational wave candidates found during veto times will be sub- 
jected to greater scrutiny than those found during other times. 



TABLE I: Summary of the amount of data analyzed in our various 
data sets. In S3 we only analyze data from the LHO detectors when 
both HI and H2 are in science mode. Around 9% of the data is 
classified as playground data and is used to tune the parameters of 
the search. 



Data type 


Total analyzed (hours) 


Non-playground (hours) 


H1-H2 


604 


548 


H1-H2-L1 


184 


167 



We can compare the sensitivities of the LIGO detectors by 
measuring the horizon distance of a particular source — this 
is the distance to which an optimally oriented source can be 
observed with SNR = 8. In Fig. |5]we plot the horizon dis- 
tance of a (2, 16)Mq binary. This choice of component mass 
reflects that the template bank used for this search (see Sec- 
tion |V]i achieves highest matches for asymmetric binaries. In 
Fig. 1 of 0] the horizon distance for a range of symmetric 
binaries is shown. 
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FIG. 5: Distance to which an optimally oriented non-spinning 
(2, 16)AfQ binary can be detected with SNR — 8 throughout S3. 
For systems with spinning components, the horizon distance would 
be equal or less than what is shown in this figure since any spin- 
induced precession would cause the system to become less than op- 
timally oriented and therefore reduce the measured amplitude of its 
emission. We see a large improvement in the sensitivity of HI during 
this science run. 



B. Coincident analysis 

To minimize the false alarm probability we demand that a 
gravitational wave signal be observed by two or more detec- 
tors with similar parameters. In order to determine whether a 
trigger measured by one particular detector should be consid- 
ered as coincident with a trigger in another detector we define 
a set of coincidence windows. In this search we demand that 
for triggers from different detectors to be considered as coin- 
cident they must satisfy the following conditions: 

\ti-t2\ < Ati+At2+ri,2, (31) 
|i/'o,i - "00,21 < AV^O,! + A?/>o,2 (32) 
IV'3,1 - V'3,2| < Al/.3,i + A?/>3,2 (33) 

where ti, V-'o.i and tp^^i are the measured time of coalescence 
and phenomenological mass parameters measured using our 
template bank in detector i; At;, AV'o.i and A-03 ^ are our 
coincidence windows in detector i and Tij is the light travel 
time between detector locations i and j. The light travel time 
between LHO and LLO is ^ 10 ms. 

We tune our coincidence windows on the playground data 
in order to recover as many of our simulated signals as possi- 
ble whilst trying to minimize the false alarm rate. The use of 
playground data allows us to tune our search parameters with- 
out biasing the results of our full analysis. The tuning method 
used for this and the non-spinning search on S3/S4 data is 
described fully in ifsoll . Using this tuning method we find our 
coincidence windows to be equal for each detector with values 
At = 100 ms, AV'o = 40, 000 Hz'^/^ and AV'3 = 600 Hz^/^ 
The value of At used in this search is four times larger than 
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the 25 ms value used in the S3 search for non-spinning binary 
black holes 10] indicating that the estimation of arrival time of 
a gravitational waveform is less well determined in this search 
than in the non-spinning search. 

C. Combined SNR 

For coincident triggers we use a combined signal-to-noise 
ratio pc statistic based upon the individual signal-to-noise ra- 
tios Pi measured by each detector; 




In practice the parameters a and b are tuned so that the con- 
tours of false alarm generated using Eq. (|34] | separate trig- 
gers generated by software injection of simulated signals and 
background triggers as cleanly as possible HO] (see the next 
subsection for details of how we estimate the background). In 
this search we used values a = = 3 for all detectors. For 
coincident triggers found in all three detectors we use: 



D. Background Estimation 

We estimate the rate of accidental coincidences, otherwise 
known as the background or false alarm rate, for this search 
through analysis of time-shifted data. We time-shift the trig- 
gers obtained from each detector relative to each other and 
then repeat our analysis, searching for triggers that occur in 
coincidence between 2 or more of the detectors. By choos- 
ing our time-shifts to be suitably large 10 ms light travel 
time between LHO and LLO) we ensure that none of the co- 
incident triggers identified in our time-shift analysis could be 
caused by a true gravitational wave signal and can therefore 
be used as an estimate of the rate of accidental coincidences. 
In practice we leave HI data unshifted and time-shift H2 and 
LI by increments of 10 and 5 s respectively. In this search, we 
analyzed 100 sets of time-shifted data (50 forward shifts and 
50 backward shifts). For clarity we will use the term in-time 
to mean triggers which have not been time-shifted. 

VII. VETOES 

A. Instrument based vetoes 

We are able to veto some background triggers by observing 
correlation between the gravitational wave channel (AS_Q) of 
a particular detector and one or more of its auxiliary chan- 
nels which monitor the local physical environment. Since we 
would not expect a true gravitational wave signal to excite the 
auxiliary channels, we will treat as suspicious any excitation 
in the gravitational wave channel that is coincident in time 



with excitations in the auxiliary channels. A list of auxiliary 
channels found to effectively veto spurious (non-gravitational 
wave coincident triggers) were identified and used for all S3 
searches ifsill . Additional vetoes based upon other auxiliary 
channels were considered but were subsequently abandoned 
because the total amount of data these channels would have 
discounted, known as the dead-time, was unacceptably large. 



B. Signal based vetoes 

We can use the fact that the Hanford detectors are co- 
located to veto coincident triggers whose measured amplitude 
is not consistent between HI and H2. We check for consis- 
tency between the SNR values measured using HI and H2 
data for triggers found in coincidence. Since HI is the more 
sensitive instrument we simply required that the SNR mea- 
sured in HI be greater than that measured in H2 for an event to 
survive this veto. Since data from HI and H2 was only anal- 
ysed when both were in science mode, this veto means that 
there will be no H2-L1 coincident triggers since this would 
indicate that H2 had detected a trigger which HI was unable 
to detect. 

The veto used for the primordial black hole and binary 
neutron star searches |3l has not not been investigated for use 
in searches using detection template families (i.e., this search 
and the S2-S4 searches for non-spinning binary black holes 
SB). 



VIII. RESULTS 

In the search of the S3 LIGO data described in this pa- 
per, no triple-coincident event candidates (exceeding our pre- 
determined SNR threshold and satisfying the coincidence re- 
quirements described in Sec. IVIBt were found in triple-time 
(H1-H2-L1) data. Many double-coincident event candidates 
were found in both triple-time and double-time (H1-H2) data. 

A cumulative histogram of combined SNR for in-time and 
background coincident triggers is shown in Fig. |6] We see 
that, at the SNR threshold (i.e., the leftmost points on this 
figure), the number of in-time double-coincident triggers is 
consistent with the number of coincident triggers yielded by 
the time-shift analysis. The small excess in the number of in- 
time H1-H2 coincident triggers at higher SNRs indicates that 
there is some correlation between the LHO detectors. The 
coincident triggers contributing to this excess have been in- 
vestigated and are not believed to be caused by gravitational 
waves. Seismic activity at the Hanford site has been recorded 
throughout S3 and can cause data to become noisy simulta- 
neously in HI and H2. Coincident triggers caused by seismic 
noise will predominantly cause only in-time coincidences (al- 
though time-shift coincidences caused by two seismic events 
separated in time but shifted together can occur) leading to 
an excess of in-time coincident triggers as we have observed 
in Fig.|6] As mentioned previously, there were no coincident 
triggers observed by all three detectors. 
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A scatter plot of the SNRs measured for coincident triggers 
in H1-H2 times is shown in Fig. |7] The distribution of our 
in-time triggers is consistent with our estimation of the back- 
ground. This is also true for the double-coincident triggers 
measured in H1-H2-L1 times. 

The loudest in-time coincident trigger was observed in Hl- 
H2 when only the Hanford detectors were in science mode. 
This event candidate is measured to have SNRs of 119.3 in 
HI, 20.4 in H2 and a combined SNR of 58.3. The loudest co- 
incident triggers are subjected to systematic follow-up inves- 
tigations in which a variety of information (e.g., data quality 
at time of triggers, correlation between the detector's auxil- 
iary channels and the gravitational wave channel) is used to 
assess whether the coincident triggers could be confidently 
claimed as detection of gravitational wave events. This event 
is found at a time flagged for "conditional" vetoing. This 
means that during these times some of the detectors auxil- 
iary channels exhibited correlation with the gravitational wave 
channel (AS_Q ) and that we should be careful in how we 
treat event candidates found in these times. For this partic- 
ular coincident trigger an auxiliary channel indicated an in- 
creased numbers of dust particl es p assing through the dark 
port beam of the interferometer Iplll . Upon further investiga- 
tion it was found that this coincident trigger occurred during 
a period of seismic activity at the Hanford site and we subse- 
quently discounted this candidate as a potential gravitational 
wave event. Time-frequency images of the gravitational wave 
channel around the time of this candidate were inconsistent 
with expectations of what an inspiral signal should look like, 
further reducing the plausibility of this candidate being a true 
gravitational wave event. It is interesting, but unsurprising, 
to note that during the search for non-spinning binary black 
holes that also used S3 LIGO data, high-SNR triggers associ- 
ated with this seismic activity were also detected |7j]. Further- 
more, the 20 next loudest event candidates were also investi- 
gated and none were found to be plausible gravitational wave 
event candidates. Work is in progress to automate the follow- 
up investigative procedure and to include new techniques in- 
cluding null-stream and Markov chain Monte Carlo analysis 
for assessing the plausibility of coincident triggers as gravita- 
tional wave events. 



IX. UPPER LIMITS 
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FIG. 6: Cumulative histograms of the combined SNR, pc for in-time 
coincident triggers (triangles) and our background (crosses with one- 
sigma deviation shown) for all H1-H2 and H1-H2-L1 times within 
S3. We see a small excess in the number of in-time coincident trig- 
gers with combined SNR ~ 45. This excess was investigated and 
was caused by an excess of H1-H2 coincident triggers. Since HI 
and H2 are co-located, both detectors are affected by the same local 
disturbances (e.g., seismic activity) which contributes to the number 
of in-time coincidences but which is under-represented in time-shift 
estimates of the background. 
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Given the absence of plausible detection candidates within 
the search described above, we have calculated an upper limit 
on the rate of spinning compact object coalescence in the uni- 
verse. We quote the upper limit rate in units of yr^^Lj^Q^ 
where Lio = lO^*' Lq.b is 10^" times the blue light lumi- 
nosity of the sun. 

The absorption-corrected blue light luminosity of a galaxy 
infers its massive star formation rate which we assume scales 
with the rate of compact binary coalescence within it lf52ll . 
This assumption is well justified when the galaxies reached by 
the detector are dominated by spiral galaxies with ongoing star 
formation (e^, the Milky Way). Results papers reporting on 
SI and S2 ifisllssi Isill have quoted the upper limit in units of 



FIG. 7: Scatter plot of SNR for coincident triggers in H1-H2 times. 
The black circles represent in-time coincident triggers and the light 
colored (red) pluses represent time-shift coincident triggers that we 
use to estimate the background. Note that due to our signal based 
veto on HI/H2 SNR we see no coincident triggers with pni < Ph2. 



Milky Way Equivalent Galaxy (MWEG) which is equivalent 
to about 1.7 Lio. Upper limits on the rate of coalescences 
calculated during other searches using S3 and S4 LIGO are 
given in units of Lio 101 ■ 

The u ppe r limit calculations are based on the loudest event 
statistic i55l l56ll . which uses both the detection efficiency at 
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FIG. 8: Upper limits on the spinning binary coalescence rate per Lio 
as a function of the total mass of the binary. For this calculation, 
we have evaluated the efficiency of the search using a population of 
binary systems with mi = 1.35Mq and m2 uniformly distributed 
between 2 and 20Mq . The darker area on the plot shows the region 
excluded after marginalization over the estimated systematic errors 
whereas the lighter region shows the region excluded if these sys- 
tematic errors are ignored. The effect of marginalization is typically 
small (< 1%). The initial decrease in the upper limit corresponds to 
the increasing amplitude of the signals as total mass increases. The 
subsequent increase in upper limit is due to the counter effect that 
as total mass increases the signals become shorter and have fewer 
cycles in LlGO's frequency band of good sensitivity. 



the combined SNR of the loudest event candidate and the as- 
sociated background probability. The in-time non-playground 
data set (which we use to set the upper limit) is blinded in the 
sense that all analysis parameters are tuned (as described in 
Sees. IVIl i prior to its analysis. 

The Bayesian upper limit at a confidence level a, assuming 
a uniform prior on the rate R, is given by ll56ll 



1 



A 



1 + A 



RTCl{Pc 



(36) 

where Cl(pc, max) is the cumulative blue light luminosity to 
which we are sensitive at a given value of combined SNR 
Pc.max, T is the observation time, and A is a measure of the 
likelihood that the loudest event is consistent with being a 
signal and inconsistent with background (as estimated using 
time-shifts). We evaluate the cumulative luminosity Cl at the 
combined SNR of the loudest coincident trigger seen in this 
search, Pcmax — 58.3 (see Sec. IVIIII for discussion of this 
coincident trigger). The expression for A is 



A 



(Pc,max) | 
(Pc,max) 



(Pc,max) 



(37) 



where the derivatives are with respect to p^. Pb (p) is the prob- 
ability that all background coincident triggers (as estimated 
using time-shifts) have a combined SNR less than p. For the 



loudest event candidate in this search we find Pb = 0.23 and 
A = 0.05. In the case where the loudest event candidate is 
most Ukely due to the background A and the upper limit 
becomes 



-^90% — 



2.3 



,n 



(38) 



In the limit of zero background, i.e., the event is definitely not 
background, A — *■ oo and the numerator in Eq. ( [38] l becomes 
3.9. The observation time T is taken from TableU where we 
use the analyzed time not in the playground. This is consistent 
with our blind analysis strategy. 

In searches for systems consisting of non-spinning bodies 
efficiency is typically found as a function of its effective dis- 
tance and chirp mass |57']. For a system consisting of non- 
spinning bodies effective distance can be calculated using the 
distance to the source, its inclination with respect to the detec- 
tor and the detector's antenna response functions (see Eq. (2) 
of lH and il). For a system consisting of spinning bodies, 
its inclination with respect to a detector will evolve during the 
course of the inspiral making the calculation of effective dis- 
tance complicated. Instead, in this search we find efficiency 
and predicted source luminosity as a function of the inverse of 
the expected SNR of a source. The expected SNR is defined as 
the SNR that would be obtained for a given simulated source 
assuming we use a template that perfectly matches the emit- 
ted gravitational waveform and a detector whose noise power 
spectrum we can estimate accurately. By taking the inverse of 
the expected SNR we obtain a quantity which behaves simi- 
larly to the effective distance by taking larger values for sig- 
nals which are nearer and/or optimally oriented to the detector 
and thus more easily detectable and by taking smaller values 
as the signals become less detectable. 

Following the tests of the template bank (Sec. IV Cl l we 
also know that the efficiency at which we are able to detect 
sources will depend on their spins as well as their effective dis- 
tance and component masses. In this upper limit calculation 
we assess the efficiency of the search using software injec- 
tion of simulated signals representing a population of sources 
with spins randomized so that i) the spin magnitude of each 
of the compact objects is distributed uniformly in the range 
< X < 1 and ii) the direction of compact object's spin is uni- 
formly distributed on the surface of a sphere. The distances of 
the simulated sources are chosen uniformly on a logarithmic 
scale. The sky-positions and initial polarization and inclina- 
tion angles of the simulated sources are all chosen randomly 
and to be uniformly distributed on the surface of a sphere. 
We evaluated the efficiency of this search for masses in the 
range 1.0 Mq < mi,TO2 < 20.0 Mq. During S3, LIGO's 
efficiency was dominated by sources within the Milky Way 
for which detection efficiency was high across the entire mass 
range investigated due to the proximity of these sources to 
Earth. We also had some detection efficiency for binaries in 
M31 and M33. 

The cumulative luminosity Cl{Pc) can be obtained by gen- 
erating a population of simulated signals using information on 
the observed distribution of sources from standard astronomy 
catalogs. We use a model based on ifssll for the distribution of 
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blue light luminosity throughout the nearby Universe. We use 
software injection of simulated signals (the target waveforms 
described in SecHllli to evaluate the efficiency Z for observing 
an event with combined SNR greater than p^, as a function of 
the source's expected SNR. We then integrate E times the pre- 
dicted source luminosity L as a function of expected SNR and 



where Dp is the distance measure equal to the inverse of the 
expected SNR, at LHO (H) or LLO (L). As mentioned earlier, 
we evaluate Cl at Pcmax = 58.3. The cumulative luminosity 
was measured to be ~ 1.9Lio and is dominated by the Milky 
Way (1.7Lio) with the remainder made up by M31 and M33. 

We calculate the upper limit on the rate of coalescence for 
proto-typical NS-BH binaries with masses m\ ~ 1.35 Af© 
and TO2 ^ 5 A/q. These values correspond to a population 
of NS-BH binaries with component masses similar to those 
used to assess the NS-NS and BH-BH upper limits in |01. 
To calculate this upper limit we evaluate the efficiency of our 
search using binaries with a Gaussian mass distribution with 
means, m\ = 1.35 Mq and m-i = 5 A/q with standard de- 
viations (Ti = 0.04 Mq and a-i = \ Mq. These efficiencies 
are measured with simulated injected signals, using the same 
pipeline we used to find our candidates, counting the number 
of injections detected with SNR above Pcmax, and the num- 
ber missed. Assuming a Gaussian distribution of masses, we 
obtain an upper limit of 7?.go% = 15.8 yr~^ Lio^^. The up- 
per limit calculation takes into account the possible systematic 
uncertainties which arise in this search, which are described 
in some detail in iIstIi . and we will follow the analysis pre- 
sented there to calculate the systematic errors for the above 
result. The most significant effects are due to the possible cal- 
ibration inaccuracies of the detectors (estimated using hard- 
ware injections of simulated signals) and the finite number of 
Monte Carlo injections performed. 

We must also evaluate the systematic errors associated with 
the chosen astrophysical model of potential sources within the 
galaxy. We obtain upper limits on the rate after marginal- 
ization over the estimated systematic errors, as described 
in lf56[|57ll . After marginalization over these errors we obtain 
an upper limit of 7^90% = 15.9 yr^^ Lio^^. 

We also calculate upper limits for a range of binary sys- 
tems with mi = 1.35Af0 and m2 uniformly distributed be- 
tween 2 and 2OAf0. These upper limits, both before and after 
marginalization are shown in Fig.|8] 



X. CONCLUSIONS 

In this paper we have described the first search for gravita- 
tional waves emitted during the inspiral of compact binaries 



mass. Since a binary system will generally have slightly dif- 
ferent orientations with respect to the two LIGO observatory 
sites, the detectors at the two sites will both measure slightly 
different expected SNRs. The source's luminosity and the effi- 
ciency with which it is detected are functions of both expected 
SNRs, and the integration needed is two-dimensional: 



(39) 

I 

with spinning component bodies, which was carried out us- 
ing data taken during the third LIGO science run. Interaction 
between the binary's orbital angular momentum and the spin 
angular momenta of its components will cause precession of 
its orbital plane resulting in the modulation of the observed 
gravitational wave. 

This search uses a detection template family designed spe- 
cially to capture the spin-induced modulations of the gravi- 
tational waveform which could have resulted in them being 
missed by other searches targeted at non-spinning systems. 
The search pipeline used to carry out this and the other recent 
inspiral searches has been significantly irnproved since S2 and 
is fully described in a companion paper \]^. 

There were no plausible gravitational wave event candi- 
dates detected within the 788 hours of S3 data analyzed. The 
upper limit on the rate of coalescence for prototypical NS-BH 
binaries with spinning component bodies was calculated to be 
^90% = 15.9 yr^^ Lio^^ once errors had been marginalized 
over 

The S5 LIGO data is greatly improved in sensitivity and 
observation time than previous data sets. There is consider- 
able work in progress to further improve the techniques used 
to search for binaries with spinning component bodies us- 
ing the S5 data. This includes development of an improved 
parameter-space metric which does not depend on the strong 
modulation approximation and allows us to search a larger 
region of parameter space, including binaries with compara- 
ble masses. Preparation for another search for binaries with 
spinning component bodies using a template family described 
by physical (rather than phenomenological) parameters ll44ll 
is also underway. 
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